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We show that the magnetic anisotropy energy of antiferromagnetic ferrihydrite depends on the 
square root of the nanoparticles volume, using a method based on the analysis of statistical distribu- 
tions. The size distribution was obtained by transmission electron microscopy, and the anisotropy 
energy distributions were obtained from ac magnetic susceptibility and magnetic relaxation. The 
square root dependence corresponds to random local anisotropy, whose average is given by its vari- 
ance, and can be understood in terms of the recently proposed single phase homogeneous structure 
of ferrihydrite. 



I. INTRODUCTION 

The relation between structural and magnetic proper- 
ties is of importance from the point of view of applied 
and fundamental research. This relation is not straight- 
forward in systems with antiferromagnetic (AF) interac- 
tions, reduced dimensionality or size such as nanoparti- 
cles. In these systems, surface effects and disorder play an 
important role and therefore deviations to the superpara- 
magnetic (SP) canonic behavior are expected. Such ef- 
fects change the relation between anisotropy energy, Ea, 
volume, V , and magnetic moment, /x, found for typical 
SP systems for which Ea and /i are proportional to V . 
This is also the case of ultrathin films, where anisotropy 
energy is proportional to surface area, leading to per- 
pendicular magnetization. Important contribution from 
surface anisotropy is also found in SP nanoparticles with 
ferromagnetic interactions, where surface atoms consti- 
tute a relevant fraction of the total atoms [l| . Another 
example of non-proportionality between Ea and V is two- 
dimensional Co nanostructures, where Ea was found to 
depend on the perimeter [2|. 

The deviations to the proportionality between V and /i 
found in AF nanoparticles are associated to the fact that, 
in these systems, the net magnetic moment arises from 
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the uncompensated and/or canted moments, ^um that 
can be present at the surface, throughout the volume, or 
both. The relation between /i„„ and V reflect the origin 
of the moments. In particular, Neel has shown that /i.u„ is 
proportional to with q ~ 1/2 for moments randomly 
distributed in the volume, 1/3 for moments randomly dis- 
tributed in the surface and 2/3 for moments distributed 
throughout the surface in active planes [3]. In ferritin, a 
protein where Fe^+ is stored as ferrihydrite, q was esti- 
mated to be of the order of g = 1/2 or to be between 
1/2 and 1/3 [1 based on magnetization measurements. 
These values were obtained either by using a system with 
a given size and estimating the power relation between 
the total number of ions and the equivalent uncompen- 
sated number of ions [3j j or by the usual comparison 
of systems with different average sizes [sj . The latter ap- 
proach is limited by the possibility of synthesizing identi- 
cal systems with different average volumes, that usually 
covers less that one order of magnitude. An alternative 
approach that takes advantage from the size distribution 
is developed here. A sample with a wide distribution can 
be regarded as one system containing a set of different 
average sizes. An approach reminiscent of this approach 
was first used by Luis and co-workers for the determi- 
nation of the origin of magnetic anisotropy in gaussian 
size-distributed Co nanoparticles [l|. They concluded 
that surface anisotropy has an important contribution, 
since the Ea distribution is narrower than the V distri- 
bution. The effect of size distributions on the magnetic 
properties was later used to study two-dimensional Co 
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structures by Rusponi et al. [g]. The idea was based on 
the fact that the shape of the in-phase component of the 
ac susceptibihty x' was critically dependent on the cho- 
sen distribution, namely surface, perimeter and perime- 
ter plus surface distributions. The authors concluded 
that perimeter atoms were those relevant to the reversal 
process in the Co structures, i. e., Ea depends on the 
perimeter. Gilles and co-workers have also tried to use 
susceptibility curves to obtain the relation between 
and V in ferritin Q but found that their experimental 
curves were not very sensitive to the particular shape of 
distribution nor the value of q [7]. In a different context, 
the luminosity and the size distributions of rare-earth- 
doped nanoparticles were used to establish the relation 
between luminosity and size through the size dependent 
optical detection probability [1] . 

In this report we show that lognormal distributed 
nanoparticle samples are particularly useful to study the 
relation between a physical property and size. This is 
based in the fact that when two physical quantities are 
related by a power function, the power factor can be 
readily obtained by comparing the respective lognormal 
deviations, due to reproductive properties of the lognor- 
mal distribution function Q . Moreover we generalize this 
concept to any distribution function. This method is of 
general use and can be simply applied in cases where the 
size of the system determined a given physical property 
by a power law relation, as in the optical properties of 
quantum dots To*!. In the present context of magnetism, 
we apply this approach to AF ferrihydrite nanoparticles 
grown in a hybrid matrix to investigate the relation be- 
tween Ea and V. 



II. MODEL 

A. Relation between distributed quantities 

As pointed out in the previous section, in AF nanopar- 
ticles there is no a priori established relation between V 
and fiun- At the same time, Ea and V can also be not 
proportional. One may however expect that, in general 



Ea = aVP 



(1) 



where p can be different than 1. In a given situation 
where the average values {V) and (Ea) of one sample are 
known it is impossible to determine a and p simultane- 
ously. Their determination is usually carried out compar- 
ing samples with different {V) , considering that a and p 
are constant in all samples. Here we show how to de- 
termine a and p using magnetic studies on one lognor- 
mal distributed sample. The probability distribution of 
Ea, g{Ea), is a function of the V probability distribution 
9{V): 



If fiy) is a lognormal distribution function with param- 
eters sv and ny defined as: 



f{V) = 



1 



Vsv^/2^ 
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with: 



He = 

SE 



psv 



(5) 



This means that if V presents a lognormal distribution, 
all other physical quantities that can be related to V by 
a power relation (Eq. [1} are also lognormal distributed. 
More important, when comparing V and Ea, the ratio 
between the distribution parameters s is a direct mea- 
sure of the power p, while the relationship between n 
values gives information about a. Therefore, the relation 
between V and Ea in one sample can be quantitatively 
derived knowing the lognormal distribution of V and Ea- 
As one might expect, this method can be used to de- 
termine the relation between any two physical quantities 
related by a precise power law similar to Eq. [TJ 

The relations expressed in Eqs. [4] and [5] are a par- 
ticular case of the reproductive properties of the lognor- 
mal distribution function . In general, if Xi are inde- 
pendent random variables having lognormal distribution 
functions with parameters ni and St (as defined in Eq. [3]) , 
their product ^ = c H (with bi and c > being con- 
stants) is also lognormal distributed, with sy — ^ hiSi 
and ny — cJ2''T-i' [9]- I^i general, reproductive prop- 
erties can be used in the analysis of an output whose 
inputs are lognormal distributed, as for instance in quan- 
titative analysis of human information processing during 
psychophysical tasks ll|. However, to the best of our 



knowledge, here is the first time that they are used in 
the context of physical properties of nanoparticles. 

Although many physical properties of interest as size 
are often lognormal distributed many others are better 
characterized by other functions. This is the case of the 
anisotropy ener gy, which is often described by a gamma 
distribution [T^ Il3l. [l^ [Tsj . The gamma function can be 
expressed by: 



Via) 



■ exp 



(6) 



g{Ea) = f[V)/{dEa/dV) 



(2) 



with the average of x given by ah and the variance by 
tj = ab"^. For a > 1, the gamma distribution is similar 
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to the lognormal function, so that the use of the lat- 
ter function in the case where the gamma distribution is 
more suitable may be a good approximation. Therefore 
the use of Eq. may also be a good approximation to 
find a and p. These values may also be found in the 
general case of a different or an unknown distribution, 
by searching for a scaling plot or numerically [? ] but, 
as seen, the validity of the lognormal distribution makes 
this task quite straightforward. 



B. Anisotropy energy distribution from ac 
susceptibility and viscosity measurements 

The out-of-phase component of the ac susceptibility 
x" is usually used to obtain the anisotropy energy bar- 
rier distribution of different nanoparticles systems and is 
given by 0, El, Q E Ofi] : 

x"(/,r) !^fc5Tln(l/(/ro))/(K) (7) 

where tq is a microscopic characteristic time and Ea = 
fcBrin(l/(/ro)) is the activation energy of the particles 
having r equal to the characteristic time of measurement 
1//. It follows from Eq. [7]that x" is a function of Ea and 
that therefore curves taken at different frequencies should 
scale when plotted against E^. At the same time, ^' jT 
is a measure of the anisotropy energy distribution f{Ec). 
In Eq. [7] it is considered that the particles contributing 
to x" ci give n / and T are mainly those with energy 
equal to Ea [I4I and that the parallel susceptibility is well 
approximated by the equilibrium susceptibility [l2l (i- e. 
1/ ln(l/(/To)) ^ 1). It is also considered that dipolar 
interactions arc negligible. 

Measurements of the magnetization as a function of 
time t (viscosity measurements) at temperatures below 
the blocking temperature Tg are a complementary way 
to investigate the anisotropy energy barrier distribution 
of different nanoparticles systems fl|, [3 El) US] , includ- 
ing ferritin [2l|. With such measurements it is possible 
to determine the magnetic viscosity, 5", defined as the 
change in magnetization with hi{t) of a system held un- 
der a constant applied magnetic field, h, and may be 
written as: 

S{t, T) = ^ = kBTM,,f{Ea) (8) 

considering that the function (t/r) exp(— ^/r) is narrower 
than the distribution function f{Ea) [l8|. In ferromag- 
netic materials M^q/h = Mg/3K, where Ms is the satu- 
ration magnetization and K the anisotropy constant. It 
follows directly from Eq[8]that S/T is proportional to the 
anisotropy energy distribution, f{Ea), in analogy with 
x" /T. In fact, S and x" are probing the same energy 
barrier at different time scales. 



III. EXPERIMENTAL DETAILS 

Ferrihydrite is a low-crystalline AF iron oxide- 
hydroxide that typically forms after ra pid hydrolysis of 
iron at low pH and low temperatures [22| . The struc- 
ture of ferrihydrite with domain sizes ranging from 2 to 
6 nm was recently described as a single phase, based on 
the packing of clusters, constituted by one tetrahedrally 
coordinated Fe atom surrounded by 12 octahedraUy co- 
ordinated Fe atoms [1^. The cell dimensions and site 
occupancies change slightly and systematically with av- 
erage domain size, reflecting some disorder and relaxation 
effects. This picture extends homogeneously to the sur- 
face of the domains. This model contrasts with previous 
ones, where multiple structural phases were considered 
[2^ l25j , and the existence of tetrahedrally coordinated 
Fe atoms was a matter of debate [1^ [131 . 

The synthesis of the ferrihydrite nanoparticles in the 
organic-inorganic matrix (termed di-ureasil) has been de- 
scribed elsewhere [11]. The particles are precipitated by 
thermal treatment at 80 °C, after the incorporation of 
iron nitrate in the matrix. The sample studied here has 
an iron concentration of 2.1 wt% and was structurally 
characterized in detail in Ref. [29|. Mossbauer spec- 
troscopy was measured at selected temperatures between 
4.2 K and 40 K. A conventional constant-acceleration 
spectrometer was used in transmission geometry with a 
^''Co/Rh source, using a a-Fe foil at room temperature to 
calibrate isomer shifts and velocity scale. Ac and dc mag- 
netic measurements were performed in a Quantum De- 
sign superconducting quantum interference device mag- 
netometer. 



IV. RESULTS AND DISCUSSION 

A. Relationship between anisotropy energy and 
size 

The Fourier transform high resolution transmission 
electron microscopy images (FT-HRTEM) and XRD 
diffraction patterns show the existence of low crystalline 
6-linc ferrihydrite nanoparticles. The nanoparticles are 
homogeneously distributed, separated from each other, 
and have globular habit. The size (diameter, D) distri- 
bution can be described by a lognormal function, with 
no = 4.7 ± 0.2 nm and deviation sd = 0.43 ± 0.05 [H 
(see Fig. [5]). As expected from the reproductive proper- 
ties, a lognormal size distribution results in a lognormal 
volume distribution. 

The in-phase ac susceptibility, x', is frequency inde- 
pendent above Tp = 30 K. The maxima of x' follow a 
Neel-Arrhenius relation: 

The extrapolated Tq is of the order of 10~^^ s, as found 
in non-interacting/very weakly interacting nanoparticles 
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|3(j |. As dipolar interactions become relevant, the ex- 
trapolated tq increases. For instance, similar ferrihy- 
drite/hybrid matrix composites with more concentrated 
fcrrihydrite nanoparticles (6.5% of iron in weight), and 
thus relevant dipolar interactions, have extrapolated tq « 
10-" s. 

Another evidence of the existence of negligible dipolar 
interactions is given by Mossbauer spectroscopy results, 
since interacting systems have a collapsed V-shaped pat- 
tern [Sll,!!!,!!!]. For temperatures around the spectra 
can be described by the simple sum of a sextet distribu- 
tion and a doublet and no signs of a collapsed magnetic 
hyperfine field pattern. On the other hand, such collapse 
is observed in the ferrihydrite/hybrid matrix sample with 
6.5% of iron, where dipolar interactions are expected to 
be relevant. At 4.2 K, the Mossbauer spectrum of the 
sample here studied (2.1% of iron) shows a sextet, with 
a hyperfine field B^f — 48 T. This is characteristic of 
fcrrihydrite nanoparticles low crystallinity, in accordance 
with the FT-HRTEM and XRD results. 

As described in Sec. HTBI x" /T and S/T constitute 
a direct measure of the anisotropy energy distribution, 
observed at different time scales. In Fig. [T] we can 
observe that the distribution obtained from x" /T and 
S/T fairly superimpose, meaning that Eq. [7] and [8] are 
good approximations. Both x" /T and S/T curves are 
well fitted by a gamma distribution function, with a=3.3 
and b=53 (Fig. [T]). As expected for a > 1, the curves 
can also be satisfactorily fitted to a lognormal function, 
with s^/T = 0.61 ± 0.02 and {x/T) = 170 ± 4 K, and 
Ss/T = 0.65 ± 0.02 and {S/T) = 176 ± 4 K, respectively. 
We therefore consider se = 0.63 ± 0.04 from the average 
of s-y^/T and ss/t- Since sd = 0.43 ± 0.05, and using Eq. 
[5] we directly obtain the power relation between Ea and 
D,p' = 1.5±0.2, which corresponds top' = 3/2, so that: 
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FIG. 1: Anisotropy energy distributions obtained by the out- 
of-phase component of ac susceptibility {x" /T) and viscosity 
(S/T) measurements. Lognormal and gamma distribution fits 
to x" /T data are shown. 
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Eq. [5] can be further used to determine the propor- 
tionahty between Ea and D^/"^, a' = 18 Knm~^/^. As 
expected from the above equation, we observe that, in 
a {Ea/a'f^ scale, both distributions of x" /T and S/T 
superimpose to the diameter distribution (Fig. [2|). This 
is a confirmation that describing x/T and S'/T by a log- 
normal function is a good approximation for the identi- 
fication of p and a. Eq. [10] can be rewritten in terms of 
the particle volume as: 



(11) 



This means that the anisotropy barriers are randomly 
distributed in volume. In each particle, the effective value 
oiEa is given by the fluctuation of local Ea- This requires 
that the local Ea is a random homogeneous quantity. 
Such homogeneity is supported by the structure model, 
since it is composed of a single phase with in-volume 
defects, where we can expected Ea to be locally different. 

From Eq. [TT]it is still possible to determine an effective 
anisotropy energy per volume, Keff, that increases with 



FIG. 2: Diameter distribution determined by TEM compared 
to the anisotropy energy distributions obtained by x" /T 
and S/T in a {Ea/a')^^^ scale, showing the scaling between 
(Ea/a)^^^ and D. Using other powers p' such as 3, 2, or 1 
gives unsatisfactory scaling. 



decreasing V, following V^^^^. In the 1 -10 nm D range, 
Keff ranges from 4.7 x 10^ to 1.5 x 10** J/m"^ which are 
of the order of those found in the literature 0, [34[ . For 
the average size of the sample, Keff = 2.9 x lO'' J/m^. 



B. Relationship between magnetic moment and 
size 

Unlike the case of Ea , there is no direct measurement 
of the fj,un distribution. A way to obtain this distribution 
is to model the dependence of the magnetization with the 
field M{H) to a given function of /i„„ considering a 
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distribution. A function usually applied to model M(H) 
of nanoparticles is the Langevin function 0, S IS]- 
This is a good approximation when surface effects and 
anisotropy are negligible. Anisotropy effects are expected 
to be relevant in AF nanoparticles due to coupling be- 
tween Hun and the AF axis 0, [lEl- In AF nanoparticles, 
anisotropy effects have been taken into account using a 
Neel (Ising-like) model, considering that can have 
only the AF axis direction [?[. On the other hand, recent 
simulations show that M{H) is greatly affected by sur- 
face effects, such that a one-spin approach as considered 
in the Langevin or Neel functions are crude approxima- 
tions (361] . 

Despite this situation, we have previously modelled 
M (H) using a Langevin distributed function '29*1 and 
found that the parameter s of the /i„„ lognormal distri- 
bution is — 0.9, so that q — 0.7±0.1 and /i„„ cx V^^^. 
We note that the value of q here derived is different to 
that estimated in Ref. [2^ comparing the average val- 
ues of the equivalent number of uncompensated ions and 
the total number of ions {q = 1/3). Both q values are 
obtained after the same fit procedure performed on the 
same M{H) curves. The only difference is the approach 
for deriving q: using average values of the uncompen- 
sated moment and size or using the information about 
the distribution of both. This is an example of how the 
use of averages may lead to inaccurate estimations, since 
the pre-factor of the power law cannot be ignored. In this 
scenario the uncompensated moments were to lie on the 
particles surface, despite the fact that the energy barriers 
associated to the uncompensated moments are randomly 
distributed in volume. 

At this point one should highlight that the Langevin 
distributed function may be a too crude description of 
M{H) to yield a good estimation of s^, so that a dif- 
ferent scenario is possible: having no reliable estimation 
of Sfj_ we discuss the situation where the uncompensated 
moments are associated (proportional) to the energy bar- 
riers, so that a oc V^^^, i. e. they are randomly 
distributed in volume. In fact, the uncompensated mo- 
ments arc those contributing to the Curie-like ac suscep- 
tibility and those experiencing the blocking phenomena 
associated with the onset of x" and S. Therefore uncom- 
pensated moments should be those relevant in determin- 



ing the relation between Ea and V. Within this frame- 
work, may be regarded as the equivalent volume 
that contains the ferromagnetic-like uncompensated mo- 
ments. Such relation between /i„„ and y^/^ was proposed 
for antiferromagnetic nanoparticles by Neel Q and is con- 
sistent with magnetization measurements performed on 
ferritin [i,[E@. 



V. CONCLUSIONS 

In this report we show that distributed samples can be 
used to investigate the relationship between the magnetic 
anisotropy barrier Ea and the nanoparticles volume V in 
a consistent manner. The relation is accessed by com- 
paring the parameters of the lognormal distribution of 
both physical quantities. Size distribution was obtained 
by a TEM study and Ea was obtained by two indepen- 
dent measurements: out-of-phase ac susceptibility and 
viscosity measurements. We have applied this method to 
a ferrihydrite nanoparticles system and found the rela- 
tion between Ea and V in an antiferromagnetic material: 
Ea OC V^^^. This shows that the magnetic anisotropy 
barriers are randomly distributed in the volume, in ac- 
cordance to recent structure studies. 
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